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\ coupled equations. We find only "closed" solutions. The behaviour of the 
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core lies completely inside the horizon. 
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I. INTRODUCTION 



A number of different topological defects [1] are thought to have been formed during the 
phase transitions in the early universe. Depending on the topology of the vacuum manifold 
A4 these are domain walls, strings, monopoles and textures corresponding to the homotopy 
groups TTolAi), 7ri(A^), 7r2(A^) and 773 (A^), respectively. Cosmic strings [1,2] have always 
gained a lot of interest since they arc thought to be important for the structure formation 
in the universe due to their huge energy per unit length (roughly 10^^^ for a string formed 
at GUT scale lO^^GeF). 

A classical field theory model which has string-like solutions is the Abelian Higgs model 
[3]. These solutions, also sometimes called "vortices", correspond to infinitely long objects. 
They have a core radius inverse proportional to the Higgs boson mass and magnetic flux 
tubes with radius proportional to the inverse of the gauge boson mass. Coupling the Abelian 
Higgs model minimally to gravity, the influence of the vortex on the geometry of space-time 
was investigated analytically [4]. It was shown that far away from the core of the string 
the space-time is Minkowski minus a wedge. It was also realised [5] that if the vacuum 
expectation value of the Higgs fleld is sufficiently large (corresponding to strings having 
formed at a phase transition with energy scale much higher than the GUT scale), then a 
different type of solution is possible. These so-called "supermassive strings" exist only on a 
finite interval of the radial coordinate and have a curvature singularity at the maximal value 
of the radial coordinate. The existence of further sohitions was investigated in a detailed 
numerical analysis [6,7]. It was found that the parameter space is indeed divided by the 
curve of maximal angular deficit 2n. If the deficit is smaller than 2tt, so-called "open" , i.e. 
infinitely extended solutions were found. One is the above mentioned cosmic string solution 
[4] which however has a "shadow" solution of Melvin-type for all values of the couphng 
constants. For deficit large than 27i, only the "supermassive", "closed" solutions exist [5]. 

The static solutions of the model without gauge field, so-called global strings, have also 
been studied [8-11]. Like all global defects, the global string has a long-range Goldstone 
field which leads to a divergent energy. Moreover, the global string is characterised by a 
logarithmically divergent deficit angle in contrast to the local string which has a constant 
deficit angle. The coupling to gravity in the case of the static global monopole leads to 
a singularity-free monopole solution [12] in the sense that while the energy is still linearly 
divergent, the solid deficit angle is now finite. For the static global string the corresponding 
singularity can not be removed by coupling the system to gravity and only the assumption 
that the metric be time-dependent removes the singularity [13]. 

Since a number of astrophysical observations like e.g. the measurement of redshifts 
of Type la supernovae [14] has led scientists to believe that we live in a universe with 
positive cosmological constant, the study of topological defects in de Sitter (dS) space seems 
interesting. But it also is of interest from another point of view, namely the dS/CFT 
correspondence [15]. This correspondence suggests a holographic duality between gravity in 
a (i-dimensional dS space and a conformal fleld theory (CFT) "living" on the boundary of 
the dS spacetime and thus being d — 1-dimensional. 

Recently, Nielsen-Olesen strings in the background of a 4-dimensional de Sitter spacetime 
dS^ have been studied [16]. However, to our knowledge, the full system of coupled matter 
and metric fleld equations has not been studied yet. One of the aims of this paper is the 
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investigation of exactly this point. 

Motivated by some recent work on a composite system of a global and local monopole 
in curved space-time [17-19], we investigate the composite system of a global and Nielsen- 
Olesen string in a curved space-time with cosmological constant as well. 

Our paper is organised as follows: we give the model and static, cylindrically symmetric 
Ansatz in Section 2. We give the equations of motion in Section 3. We discuss the pure 
Nielscn-OIcsen solutions in Section 4, the global string solutions in Section 5 and the com- 
posite system of a global and Nielsen-Olesen string in Section 6. We give our summary in 
Section 7. 

II. THE MODEL 

The model which describes a gravitating Nielsen-Olesen string interacting with a global 
one in the presence of a non- vanishing cosmological constant is given by the following action: 

S = J d^X^f^ (l6^^^ ~ '^^^ ^ ^ ^global + Anter) (l) 

where R is the Ricci scalar, G denotes Newton's constant and A is the cosmological constant. 
The Lagrangian of the Abelian Higgs model is given by [3] : 

Cno = \D,<t>{D^^<t>r - \f,.F>^^ - ^ (#* - r){f (2) 

with the covariant derivative = — ieA^ and the field strength F^y — d^A^ — d^A^ of 
the U(l) gauge potential A^ with coupling constant e. is a complex scalar field (the Higgs 
field) with vacuum expectation value rji and self-coupling constant Ai. The Lagrangian of 
the global string reads [1] : 

C,iobai = Id.xd^x* - ^ {XX* - ril) ' (3) 

where x is a complex scalar field (the Goldstonc field) with vacuum expectation value r]2 
and self-coupling A2. Finally, following [17] we introduce an extra potential which couples 
(with coupling constant A3) the two sectors of the model directly to each other: 

Ante. = (00* - r,t) (xX* - Vl) (4) 

Without this term, the global and local string would be coupled only indirectly over gravity. 
In this paper we will use units which h — c — 1. 

A. The Ansatz 

In the following we shall analyse the classical equations of motion associated with the 
above system. In order to do that, let us write down the matter and gravitational fields as 
shown below. The most general, cylindrically symmetric line element invariant under boosts 
along the direction is: 
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ds^ = N''{p)dt'' - dp"" - L\p)dif'' - N\p)dz^ . (5) 
The non-vanishing components of the Einstein tensor G^,^ then read: 

Gtt = -C.. = J {LdppN + d,Nd,L + NdppL) , 

d N LP' 
Grr = ^ (29,L7V + dpNL) , G^^^ — {2Nd,,N + {d,Nf) , (6) 

where dp denotes the derivative with respect to p. 
For the matter and gauge fields, we have: 

0(p,(^)=7yiMp)e*"^ (7) 
X(p,</:')=m/(p)e'"^^ (8) 

A^dx^ = ^(n - P{p))d^ . (9) 

n and m are integers indexing the vorticity of the Higgs and Goldstone fields, respectively, 

around the 2;— axis. 

Substituting the above configurations into the matter Lagrangian density Cm = + 
^global + Winter, we obtain: 

- ^ifiP) - - ^mp) - imp) - g^) . (lo) 



III. EQUATIONS OF MOTION 

We define the following dimensionless variable and function: 

X = \[h.r]ip , L{x) = L{p)r]i^fxl . (11) 
Then, the total Lagrangian only depends on the following dimensionless coupling constants 

7 = 87rG77i ,« = eVAi , q ^ , A ^ , ^ ^ , i = 1,2,3. (12) 

Varying (1) with respect to the matter fields and metric functions, we obtain a system 
of five non-linear differential equations. The Euler-Lagrange equations for the matter field 
functions read: 
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{N^Lh'y 



n 



hP^ + h[h' _ 1) + p^^h{f' _ q^) , 



(13) 



mi 



L2 



+ Pifif - q') + Pifih' - 1) , 



(14) 



m [ L 



ah^P , 



while the Einstein equations 

G^i. + Ag^^ = 'jT^^ , ij,,u = t,x,(p,z 



read: 



{LNNJ 



= -A + 7 



2q;L2 



2\2 



and 
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(iV^L 
iV2L 



= -A - 7 



n^h^(x)P\x) m^f(x) n\P'{x)f I,,,/ n 

+ \ + ^ \ + 7(^ (2^) - 1) 



L2(x) 



L2(x) 2aL2(a;) 4 
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+ ^(Ax) - g^)^ + f (/.^(x) - l)(Ax) - g^) 



Moreover, defining u — — N'^L we get the following equation: 



u"{x) 
u{x) 



= -3A - 7 



n^h^{x)P^{x) m^f{x) n\P'{x)f 



L-'ix) 



+ 



L^{x) 2aL'^{x) 



+ ^(/.^(a;)-l)^ + ^(f(x)-gT 
+ ^Ji^h\x)-l){f{x)-q') 



The prime now denotes the derivatives with respect to x. 



(15) 



(16) 



(17) 



(18) 



(19) 



A. Boundary Conditions 

The requirement of regularity at the origin leads to the following boundary conditions: 

hiO) = 0, /(O) = , P(0) = n (20) 

for the matter fields and 
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A^(0) = 1, A^'(O) = 0, L(0) = , L'(0) = 1 



(21) 



for the metric fields. Since a cosmological horizon appears naturally in de Sitter space, we 
integrate the equations only up to this value of the coordinate x, x — xq. In order for the 
core of the local string to lie completely within the horizon we require: 

h{x = xo) = 1, fix = xo)=q, P{x = xo) = . (22) 

Note that due to the fact that the x interval is finite, this is not (like in asymptotically 
flat space) a necessary condition for finite energy solutions. However, we have chosen these 
boundary conditions such that the energy-momentum tensor vanishes at a; = a;o- In addition, 
the limit A — > which leads to — > oo can be taken with these boundary conditions. 



IV. NIELSEN-OLESEN STRINGS IN DE SITTER SPACE 

First, we are interested in the case of the pure Nielsen-Olesen string. This corresponds 
to setting f{x) = and g = in the previous equations. 

A. Vacuum solution 

For the case of the pure gauge string, there is a vacuum solution of the equations. Setting 
P{x) — and h{x) — 1, we find from (19), that: 



N'^{x)L{x) = Asin(V3Aa;) + Scos(V3Ax) , A,B constants . (23) 
Using the boundary conditions (21), we find the following solution: 

N'^{x)L{x) = -i^sin(v^x) . (24) 
V3A 

This then can be put into (17) and (18) and we find the solutions: 

N{x) = cos2/3(y3A-) (25) 
2 

and 

02/3 

L{x) = ^[sin(V3Ax)]^/3[tan(V3A-)]2/3 (26) 
V3A 2 

where the coefficients again result from the boundary conditions (21). The first zero of N{x) 
lies at Xq = tt/VsA. At the same time, L(a; — > Xq) — > cxo. This is the cosmological horizon 
of the vacuum solution. If we expand the metric functions around this horizon, we obtain: 

Nix xl) ^ i-^m^ - xlf" + ... (27) 

and 
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L{x xl) ^ {^)-"\x - xl)-'" + .. . (28) 

This has the behaviour of a so-called Kasner solution [20] : 

ds" = {kpf'^de - dp" - C^kpf^'-'^p^d^^ - {kpfdz" (29) 

with a = c — 2/3, 6 = — 1/3, A; = C = 1. These type of "closed" solutions have 

been found previously [5] for the case A = and were called "supermassive" strings. When 
calculating the Kretschmann scalar K = R^'''^'' R^^pa- one obtains [5] that K oc {x — Xq)^'^ 
and thus the solution has indeed a curvature singularity ai x — Xq. Remarkable is that in the 
case of A = 0, these type of solutions only appear for sufficiently high vacuum expectation 
values (vev) of the Higgs field corresponding to strings having formed at energy scales much 
higher than the GUT scale [5]. For smaller values of the vev no singularity appears and the 
solutions exist on the full interval [0 : oo[. Accordingly, the numerical study showed [7] that 
these solutions exist for a 7 > 7cr. 



B. Numerical results 

Subject to the boundary conditions (21), we have studied the coupled system of equations 
numerically. 

First, we fixed a and A to study the influence of the gravitational coupling 7 on the 
solutions. We determined the value of x at which the metric function N{x) vanishes, i.e. 
N{x — Xo) — 0. Our results for a — 1.0, A = 0.005 and n = 1, 2 are shown in Fig. 1. 
As expected the value of xq decreases with increasing gravitational coupling 7. Moreover, 
we observe a steep decrease in xq for a relatively small range of 7. We have only plotted 
results for 7's corresponding to a^o > 5 since for large 7's the numerics becomes increasingly 
difficult. The reason for this is indicated in Fig. 2, where we show the profiles of the metric 
functions N(x), L{x) as well as those of the matter field functions P{x) and h{x) for a — 1.0, 
A = 0.005, n = 1 and two different choices of 7. For 7 = 1.5, the value of x at which the 
matter field functions reach their asymptotic values and 1, respectively, is much smaller 
than the value of xq. This means that the horizon clearly lies outside the core of the string. 
For 7 = 1.7, however, the situation is different. The value of x at which h{x) reaches 1 
is roughly equal to Xq, while P{x) seems to be still greater than on the plot we present. 
This is due to the fact that the "real" solution would have a slightly higher xq at which 
P{x = Xq) = 0. However, since L{x xq) 00, it is numerically impossible to reach 
the final solution. Nevertheless, the plot indicates that -like h{x)- P{x) just reaches its 
asymptotic value at x — xq. Thus the horizon lies very close to the core of the string. 

Then, we fixed a and 7 and determined Xq in dependence on A. Our results for a — 
7 = 0.5 and n = 1, 2 together with the location of the cosmological horizon of the vacuum 
solution, Xq, are given in Fig. 3. 

We clearly observe that the value of the cosmological horizon decreases with the increase 
of the cosmological constant, as expected. Moreover, an increase in the vorticity n leads to 
a decrease of Xq for the same A. For small A, xq of both the n — 1 and the n — 2 solution is 
very close to the corresponding Xq. This can again be explained by studying the behaviour 
of the matter field functions for varying A. As observed previously [16], we find that for 
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fixed a and 7 and increasing A, the value of the coordinate x at which the matter field 
functions reach their asymptotic values also increases, e.g. for a = 7 = 0.5, n = 2 we find 
that for A = 0.001 the value of where the gauge field function reaches P{x'^^^) = 0.1 is 
x°-i(A = 0.001) ^ 6.2, while for A = 0.005, we find x^-^{K = 0.005) ^ 6.45. This can be 
interpreted as representing a thicker string core due to an increased cosmological expansion. 
Thus the cosmological horizon lies closer and closer to the core of the string for increasing 
A and so only for small A the solution close to the cosmological horizon can be described 
by the vacuum solution. In Fig. 4, we show the profiles of N[x) and L{x) for A = 0.01, 
q; = 7 = 0.5 and n — 2 both for the Nielsen- Olesen string and the vacuum solution. Clearly, 
the solutions differ quite strongly. 

V. GLOBAL STRINGS IN DE SITTER SPACE 

This is the case of setting P{x) = n and h{x) = (which implies 771 = 0). The global 
string without a cosmological constant, i.e. A = has been studied extensively [8-11]. To 
have a good starting solution for the construction of dS global strings, we have reconstructed 
the global string without cosmological constant. We find that the metric function N{x) 
deviates very little from one, but that the quantity \N{x = 0) — N{x*)\ is increasing with 
X* 00. Moreover, L{x) grows approximately linearly with x. This is the behaviour found 
in [9], namely that outside the core of the global string the metric functions behave like 
A^^ = 1 — |ln(^), = a;^(l — 7ln(^)) with Xgc oc (g/52)~^ being the core of the global 
string in our rescaled coordinates. Moreover, we recover the behaviour of the Goldstone field 
function /(x — > 00) — q — 0{x~^). 

A. Numerical results 

For non-vanishing cosmological constant we find that the behaviour of N{x) and L{x) 
resembles that of the metric functions in the case of the Nielsen-Olesen string in de Sitter 
space (see previous section). Again, for all constants fixed and A varied, we find that the 
value at which N(x — Xq) — decreases with increasing A, e.g. for — 1, q — 0.1, 7 = 0.1, 
we find that xq{K = 10^^) ~ 180, while a:o(A = 10~^) ~ 57. At the same time, the value of 
X at which the function f{x) reaches its vev increases with increasing A which is due to the 
increased cosmological expansion thus leading to an extended string core. The behaviour 
of the function f{x) depends crucially on the cosmological constant which determines the 
cosmological horizon and the parameters P2 and q, which determine the radius of the string 
core. We will discTiss these features in more detail in the context of composite topological 
defects in the next section. 

VI. COMPOSITE SYSTEM OF GLOBAL AND NIELSEN-OLESEN STRING 

Unlike in the case of the "pure" Nielsen-Olesen string, a complete analytical analysis of 
the composite system seems to be impossible. However, some additional information can 
be gained by analysing the energy density per unit length. Before discussing our numerical 
results, we make some remarks on this point in the following section. 
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The energy density per unit length of the composite defect is given by [6,7]: 

^ = / \f^TQdxidx2 



(30) 



where is the determinant of the 2 + 1 dimensional metric ds'^ — N'^ (p) dt^ — dp^ — LF'{p) dip^ 
and Tq — —Cm is the 00-component of the energy-momentum tensor. In the cylindrical 
coordinates, we get: 



S^nrif dxN(x)L(x) 
Jo 



n 



-iP'{x)f + {h'{x)f + {f{x)f 



^ n'h'{x)P\x) ^ m^p{x) 



L2(x) L2(x) 
(3l{h{xf-l){nx)-q') 



aL'^{x) 

+ lmx)-ir+^{f{x)-q'^' 



PI 

2 



(31) 



From this, we can observe that the composite defect has a finite energy density There 
are two different scenarios now. If the product /?2g is large enough the core of the global 
string is small and the cores of both strings lie within the cosmological horizon xq. We 
can then assume that f{xo) = q, f (xq) = 0, P{xo) = 0, H^xq) = 1. From before, we 
know that L{x) goes to infinity close to the horizon, while N{x)L{x) remains finite. Thus 
the integrand of (31) tends to zero like f^'^ffj^- If the product /?2Q' is small, the core of 
the global string extends to outside the horizon Xq. Then the integrand of (31) becomes 



N{x)L{x) {f{x)f + 



m f (x) 



+ 2 (/^(2^) ~ 9^)^ • Since close to the horizon, the product 



N{x)L{x) tends to zero, this is finite. We have indeed confirmed numerically that this is 
the case. 



A. Numerical Results 

Assuming that the behaviour of the metric functions N{x) and L{x) persists in the 
presence of all (non-trivial) matter fields (which indeed our numerical analysis confirms), we 
can insert (25) and (26) into (14). We obtain that: 

f{x ^ xo) ~ g + C{x - xof^^ , C constant (32) 

We have construct the composite model solution by starting from the Nielsen-Olesen 
(NO) string and increasing the parameter q gradually from zero. Our numerical analysis 
confirms the assumption that the radius of the NO string core is smaller than that of the 
global string. Moreover, we find that the NO string always resides inside the horizon. As 
for the global string, we find that when the coupling constant (32 is small the function f{x) 
reaches the imposed expectation value /(xq) = q with a positive concavity, in particular the 
derivative f'\x=xa is non zero and our numerical solution clearly suggests that the behaviour 
f{x — >■ a;o) ~ q + C{x — XqY^^ is not fulfilled. This is demonstrated in Fig. 5, where we show 
/(x) of the composite defect for /3| = 1, « = 7 = = 1, g = 0.1, /?| = and A = 0.001. 
This suggests that for small values of the product /32g, the argument demonstrated above 
doesn't hold. 
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Increasing the value of (32 wc were able to produce solutions which seem to have the 
expected behaviour, i.e. /(O) = 0, f{x = Xq) = q and f'l^^r^^ = 0. This is shown in Fig. 5 
for a = 7 = /?2 = 1, g = 0.1, /?| = 0, A = 0.001 and /?| = 5, 10, 20, 50, respectively. 
Of course, the occurrence of the singularity Bit x — Xq renders the interpretation of the 
numerical results not hundred percent certain but we are rather confident that composite 
string defects which are regular inside the horizon exist for large enough (32- 

All these results are obtained for the case of jS^ = 0, i.e. the two defects interact with 
each other only indirectly over gravity. We have also attempted to construct solutions with 
/Ss 7^ 0. We find that only for large enough values of the quotient the solutions seem 

well behaved. If is roughly of the order of 10^ (for a — j — — 1, q — 0.1 and 

A = 0.001), the behaviour of the functions is very similar to that in the case of /^s = 0. For 
smaller than that, however, the functions f{x) and h{x) start to develop oscillations 
close to the cosmological horizon. The number of oscillations increases with the decrease of 
the quotient Thus, we conclude that directly interacting composite defects without 

the global string singularity only exist if the self-interaction of the global string is much 
larger than the interaction between the two defects. 

VII. CONCLUSION 

In this paper we have analysed both Nielsen- Olesen and global strings as well as the 
composite system of both defects in de Sitter space. When the matter fields are set equal 
to their vacuum expectation values (vev) in the case of the "pure" Nielsen- Olesen string, we 
were able to construct analytic solutions of the Einstein equations in terms of trigonometric 
functions. The metric tensor resembles that of so-called "super massive" strings which exist 
in asymptotically fiat space only for sufficiently high enough values of the vev of the Higgs 
field [5]. These were considered as being "unphysical" since they should have formed at 
energy scales high above the GUT scale. Since recent observations indicate that we live 
in a universe with positive cosmological constant, and since we find that the existence of 
our solutions is not restricted to values of the Higgs field's vev being large enough, these 
solutions might well be of relevance. 

Our numerical analysis suggests that the general behaviour of the vacuum metric persists 
in the presence of the matter fields. E.g. comparing the location of the cosmological horizon 
xo in dependence on the coupling constants for the "pure" Nielsen-Olesen string and that of 
the vacuum solution, Xq, we find that for small A and/or 7, Xq and Xq are nearly equal. For 
increasing A and/or 7, the difi^erence between the two increases. The reason for this is that 
the radius of the string core becomes comparable to the radius of the cosmological horizon 
and thus the assumption of a "vacuum" at the cosmological horizon is not valid any longer. 

Constructing the global string and the composite defect of Nielsen-Olesen and global 
string, we find that the existence of solutions without a singularity resulting from the global 
string itself depends crucially on the product f32(l which (in our rescaled coordinates) is in- 
verse proportional to the radius of the global string's core. If (for fixed q) P2 is too small, 
the core of the global string extends to outside the cosmological horizon and the function 
/ reaches its vev with positive concavity. For large enough values of ^2, the core of the 
global string lies inside the horizon and our numerical results seem to indicate that a global 
string/composite defect without the normal singularity of the global string exists. However, 
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note that the removal of the global string singularity which exists in asymptotically flat 
space can be achieved only by introducing a "new" singularity, the curvature singularity at 
the horizon. 
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Note added After finishing the manuscript, B. Linet has brought to our intention his 
paper "The static, cylindrically symmetric strings in general relativity with cosmological 
constant" [J. Math. Phys. 27 (1986), pp. 1817-1818], in which he discusses the vacuum 
solutions (25) and (26) "re-found" by us. 
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FIG. 1. The value of the dimcnsionless coordinate x at which a cosmological horizon appears, 
X = xq, is given for the Nielsen-Olesen string as function of 7 for n = 1 (dashed) and n = 2 (sohd) 
with a = 1.0 and A = 0.005. 
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FIG. 2. The profiles of the metric functions N{x), L{x) and the profiles of the matter field 
functions P{x) and h{x) are shown for a = 1.0, A = 0.005, n = 1 and 7 = 1.5 (dashed) and 7 = 1.7 
(solid), respectively. 
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FIG. 3. The value of the dimensionless coordinate x at which a cosmological horizon appears, 
X = xq, is given for the Nielsen-Olesen string as function of A for n = 1 and n = 2 with a = 7 = 0.5. 
For comparison, also the value Xg = n/VsA for the vacuum solution is given. 
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FIG. 4. The profiles of the metric functions A^(,t) and L{x) are shown for a = 7 = 0.5, A = 0.01 
and vorticity n = 2. We compare the Nielsen-Olesen solution (solid) with the corresponding vacuum 
solution (dashed). 
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